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Abstract 

The longitudinal (e, e') response function of 4 He is calculated precisely with 
full final state interaction. The explicit calculation of the four-body contin¬ 
uum states is avoided by the method of integral transforms. Precision tests 
of the response show the high level of accuracy. Non-relativistic nuclear dy¬ 
namics are used. The agreement with experimental data is very good over 
a large energy range for all considered momentum transfers (q = 300, 400, 
500 MeV/c). Only at higher q the theoretical response overestimates the 

experimental one beyond the quasi-elastic peak. 
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A new method for the calculation of the inelastic response of an IV-body system to an 
external probe is proposed in Ref. [tTJ . It allows an exact calculation without the knowledge 
of the iV-body scattering state. The high level of accuracy of the method has been shown 
for the longitudinal electron scattering responses of the nuclear two- and three-body systems 
[fj.[2j]. The real superiority of the approach, however, becomes evident when applied to a four- 
body system. In fact a solution of the four-body medium energy continuum state problem is 
presently out of reach, nonetheless four-body response functions can be reliably calculated 
as pointed out in the following. In this work we consider the important longitudinal electron 
scattering response function R r of ^He which is calculated for the transfer momenta g=300, 
400, and 500 MeV/c. For q = 500 MeV/c it is the first accurate calculation with the 
final state nuclear interaction fully taken into account. Our results are obtained within the 
framework of the non-relativistic nuclear dynamics and using the single-particle form of 
the electromagnetic operator. Such studies allow establishing the limits of validity of this 
conventional framework for the lightest tightly bound nucleus. Particularly interesting is 
the higher q region. For more than a decade there has been a lot of discussion for complex 
nuclei regarding this region. An accurate calculation for 4 He will help to shed some light in 
this range of q values. 

The idea of Ref. jl| is to calculate the response in an indirect way. First the Lorentz 
transform (LT) 



( 1 ) 


of the response function 


R(lo, q) = J2 |(n|0(q)|O)| 2 <5(u; - E n + E 0 ) 


( 2 ) 


n 


is calculated, where |0) is the ground state of the system, E 0 is the ground state energy, 
@(q) is the excitation operator, and a R > 0, oy ^ 0. The solution of the following equation 


(H — E 0 — cr R + zcr/)|T> = 0)0), 


( 3 ) 


leads directly to the LT: 
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Cf>(<7,g) = <*!*>■ 


(4) 


In a second step R(u, q) is obtained via the inversion of the transform. The solution of 
Eq. (3) is unique. Indeed, the homogenous equation has only the trivial solution because 
the hamiltonian H has only real eigenvalues. Since T has to fall off exponentially one can 
use similar techniques as for the solution of the ground state problem. Thus the extremely 
complicated asymptotic boundary condition of a four-body scattering state has not to be 
considered at all. 

In the past other integral transforms were proposed, namely Stieltjes J3| and Laplace 
transforms [|J5|]. The Laplace transforms of the longitudinal response were obtained with 
a realistic force for q = 300 and 400 MeV/c via a Green Function Monte-Carlo calculation 
(GFMC) Q. Also the Laplace transforms of the transverse response and the effects of two- 
body operators on the transforms in both longitudinal and transversal cases were considered 
via a GFMC 0. Good agreement with the transforms of the experimental data is found. 
There is, however, a fundamental problem in obtaining response functions themselves from 
these transforms. Unlike the LT they sample contributions over a large energy range. This 
results in big problems for the inversion [6]. Nevertheless the longitudinal R(u,q) of 4 He 
has been obtained by an inversion of the Laplace transform for q = 400 MeV/c J4|. The 
result is rather similar to ours in Fig. 4. We are not able to fully interpret this agreement 
since the statistical errors of a GFMC lead to an uncertainty in the inversion of the Laplace 
transform. Unfortunately, the inversion error is not estimated in ref. [j|], which in general 
can be sizeable [Q . On the contrary, for the LT inversion problems are much less important 
JUJU. Moreover, the numerical effort for the calculation of the LT seems to be much smaller 
than for the Laplace transform. However, a fair comparison can only be made when both 
calculations are performed for the same potential model. 

Our nuclear Hamiltonian includes central even potentials 


V(ij) = V 31 ( rij )P+(ij)p-(ij) + V 13 (ry)P"(y)P T + (y) 


( 5 ) 
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providing realistic description of the S -wave phase shifts np to the pion threshold. We 
construct the V 31 and V 13 potentials by modifying the complete NN interaction of Ref. 
||. The disregarded tensor force is effectively simulated via a dispersive correction (V —> 
V — V? ensor / const). The potentials obtained lead to almost the same phase shifts as in Ref. 
||. A full description of the potential will be published elsewhere 0. It describes the static 
properties of 4 He rather well leading to a binding energy of 31.3 MeV and an rms radius 
of 1.40 fm. Also the description of the elastic form factor is rather realistic up to its first 
minimum. The present ansatz for the potential will lead to results quite similar to those 
for more general nuclear forces. The three-nucleon studies undertaken so far testify to this 
opinion 0. Although more intensive the calculations with a completely realistic nuclear 
force are also quite feasible within our approach. 

In the following we describe the techniques we use for solving the dynamic equation (3). 
We seek for the solution in the form of an expansion over the correlated hyperspherical basis 
first used in Ref. [|10] . The expansion converges quickly in few-nucleon bound state problems 
iqjnfl. Our basis functions are of the form 


JRn(p) yj(LM (^) @S=0,T 


[f]n 


( 6 ) 


Here p is the hyperadius, p — (£1 + + ^I) 1 ^ & are the normalized Jacobi vectors, and 

denotes collectively eight hyperangular variables. The quantities y][^ m are hyperspheri¬ 
cal harmonics (HH) with hyperangular K and orbital L, M momentum quantum numbers. 
These HH are components p of irreducible representations [/] of the four-particle permuta¬ 
tion group 5(4). The spin-isospin functions 6 (see e.g. [Q) enter Eq. (|(J) with the same spin 
and isospin values S = 0, T — 0, and T = 1 as in the expansion of the right-hand side of 
Eq. (3). They belong to the conjugate representation [/] of 5(4). The square brackets mean 
coupling to the function antisymmetric with respect to permutations of both spatial and 
spin-isospin particle coordinates. Rn are the hyperradial functions, and J is the Jastrow 
correlation factor. 

The system of equations for the expansion coefficients is obtained by projecting Eq. 
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(3) onto the subset of functions (0[) with K up to some K max and N up to some N max . 
This system is split with respect to L, M and T values. Since L( 4 He) = 0 in our model 
the L quantum number coincides with the multipole order of the transition operator. The 
response, as well as Eq. (4), is independent of a q direction that can be chosen along the z 
axis. Only the M — 0 value gives a non-zero contribution in this case. The matrix elements 
are calculated with a Monte Carlo integration. 

The HH entering Eq. (§) are constructed by applying the convenient form, see e.g. 03, of 
the Young operators to the simple Zernike-Brinkman type HH. The multiplicities of various 
[/] representations at given K and L values are obtained as traces of the Young operators 
calculated in the Zernike-Brinkman basis fTj|]T4| . 


The hyperradial functions of the form [|TTI| . [Tof| R N (p ) ~ L 8 N (p/b) exp(— p/2b) are used. 
Here are Laguerre polynomials, and b is a scale parameter which is kept the same for 
all the a values considered and is chosen to enable sufficiently fast overall convergence. The 
results are rather insensitive to the b values. The rate of the hyperradial convergence in our 
case is lower than in the bound state calculations (e.g. |P~5| , P~5| ), and better R N can perhaps 
be found. 

The two-body correlation function f(r) entering the Jastrow factor is taken to be spin- 
independent and is chosen in a conventional way. At r < r$ it is a solution to the Schrodinger 
equation with the potential taken as the half-sum of the triplet and singlet NN forces. The 
r 0 point is chosen from the condition f'(r 0 ) = 0. At r > rq f(r ) = /(r 0 ). The kinetic energy 
matrix elements with the Jastrow factor are cast to a convenient form |10 . 


We calculate the LT of Ri,(u,q) with oq = 20 MeV. The quantities in Eqs. (2), (3) 
pertain to the center of mass system, and 


©(q) = 


1 + r 3 (k) G^jq 2 ) 1 - r 3 (fc) 
2 + G p E (q^ 2 


,«q(rt-Rc.m.) 


where G p ^ n are nucleon Sachs form factors. In order to reach convergence we choose a 
sufficiently large N max for the hyperradial functions ( N max = 20, 25, and 30 for q = 300, 400, 
and 500 MeV/c, respectively). The multipole transitions of the charge operator are taken 
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into account up to a maximal order L max . ^Frorri the evaluation of the various multipole 
contributions to the Coulomb sum rule we find that the following L max values lead to an 
exhaustion of the sum rule by more than 99%: L max = 4, 5, and 6 for q = 300, 400, and 
500 MeV/c, respectively. These L max values are adopted at solving Eq. (3). The maximal 
hyperangular order K max is taken equal to 7, only in case of L max = 6 the value 8 is used. 
This is sufficient to completely exhaust the various multipole strengths for q = 300 MeV/c. 
Even for q = 500 MeV/c one misses only a small fraction of the strength of the less important 
multipoles with L > 4 (see also discussion below). 

The results for the LT are shown in Fig. 1. Unlike Stieltjes and Laplace transforms it 
is already obvious directly from the LT that the response is governed by the quasi-elastic 
peak. The inversion is performed with the same sets of basis functions used in Refs. . 
Contrary to the nuclear two- and three-body systems, we cannot of course compare the 
R(lv, q) obtained from the inversion with a direct calculation of the response according to 
Eq. (2). Nonetheless it is possible to test the precision of the response function results. A 
first test is the separate inversion of all the various multipoles. It serves as a very important 
sum rule check, since for a given multipole one can compare the sum rule from the evaluation 
as ground state expectation value with that obtained from an explicit integration of the 
response. This check leads to very good results with relative errors of about 1% for most 
of the transitions (average errors: 1.1%, 1.0%, and 2.0% for q = 300, 400, and 500 MeV/c, 
respectively). Somewhat larger errors are found only for q = 500 MeV/c, where the less 
important higher multipoles ( L > 4) are slightly underestimated by about 3% - 4%. As 
mentioned above K max should be chosen somewhat larger for a complete exhaustion of the 
strength of these multipoles. Nonetheless we may say that the sum rule results show the 
good accuracy of our method. In Fig. 2 the isoscalar and isovector parts of the response 
function obtained from the separate inversion are shown for q = 500 MeV/c. One sees that 
almost all multipoles have the typical structure due to the quasi elastic peak. The only 
exception is the isoscalar Coulomb monopole which exhibits a peak close to threshold. For 
the two lower momentum transfers this CO peak is even more pronounced. For q = 300 
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MeV/c its height reaches already one third of the quasi-elastic peak height. The isovector 
strength is twice as large as the isoscalar one. 

Another very important check for the precision of the method is obtained by the inver¬ 
sion of the total LT. The resulting R(ui, q) should not differ from that obtained from the 
separate inversion discussed above. Before discussing these results we should mention that 
we encounter at low energy for q = 400 and 500 MeV/c similar inversion problems for the 
full LT as described in Ref. |jj. We solve this problem in a similar way as in Ref. flj, i.e. 
by separate inversions for the sum of isoscalar CO and Cl and for the sum of all remaining 
multipoles; nevertheless in the following it will be called total inversion. The total response 
functions resulting from separate and total inversions are shown in Fig. 3 for the three con¬ 
sidered momentum transfers. From the good agreement of the various curves it is evident 
that the inversion is very unproblematic. Differences between the two inversion methods are 
only found at lower energies, however they are quite unimportant. We consider the inversion 
of the total <L(cr, q) as the more accurate result, since we obtain a better fit to the calculated 
LT in the low-energy region. The total Coulomb sum rule is reproduced very precisely by 
the inversion of the total LT. We find relative errors of 0.2%, 0.4%, and 1.6% for q = 300, 
400, and 500 MeV/c, respectively. The reason for the somewhat larger error at q — 500 
MeV/c has been already discussed above. 

After having demonstrated the precision of the method we compare our results with 
experimental data. To this end we have to consider that the response function of Eq. (1) 
is defined for point particles. In order to compare with experiment we have to multiply 
R{u lab ,q) with the square of the proton charge form factor G p E (q 2 — ujf ab ), where = 
uj + g 2 /2M( 4 He). We take the dipole fit to G P E with the usual relativistic correction jl6|. In 
Fig. 4 we show our results in comparison with experimental data p~7| , |T8f . It is readily evident 
that for the lower q value of 300 MeV/c the agreement between theory and experiment is 
very good. The low-energy wings of the response at g=400 and 500 MeV/c are also in a very 
good agreement with experiment. In particular, the rather complicated threshold structure 
of the experimental R E at q = 400 and 500 MeV/c is described extremely well. Beyond 
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the quasi-elastic peak the theoretical result overestimates the experimental one somewhat 
at q — 400 MeV/c and in a more pronounced way at q — 500 MeV/c. If the experimental 
results are correct the theoretical formulation should include subnuclear and/or relativistic 
effects in order to remove the discrepancy. 

In conclusion we may say that we have successfully applied the method of Ref. [If to 
a four-body system response to an external probe with full final state interaction. This 
enabled us to calculate the accurate longitudinal response function of 4 He. We have shown 
that the results are very precise. We obtain an excellent agreement with experiment at the 
momentum transfer of 300 MeV/c as well as for the low-energy wings at q = 400 and 500 
MeV/c. At the latter q values the theoretical results overestimate the experimental ones 
beyond the quasi-elastic peak. Though somewhat more complicated a calculation with a 
fully realistic potential model can also be carried out in a similar way. The calculation of 
the transverse response with the present potential model is in progress ||. 

The authors thank INFN for having provided a dedicated work station (SUN SPARC-20) 
for the numerical calculations. One of us (V.D.E.) thanks INFN for the financial support 
over the period during which this work was carried out. 
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FIGURES 


FIG. 1. LT at q =300 (a), 400 (b), and 500 MeV/c (c) 


FIG. 2. Separate inversions of the various isoscalar (a) and isovector (b) multipoles of the LT 
(g=500 MeV/c). The various curves correspond to successive addition of multipole contributions 
from CO to C6 


FIG. 3. Response functions from total (full curves) and separate inversions (dotted curves) 

FIG. 4. Response functions from total inversions with inclusion of proton charge form factor 
(see text) in comparison to experimental data 
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